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Abstract 



We consider a new fractional order chaotic system displaying an interesting behavior. 
A necessary condition for the system to remain chaotic is derived. It is found that chaos 
exists in the system with order less than three. Using the Routh-Hurwitz and the Matignon 
stability criteria, we analyze the novel chaotic fractional order system and propose a control 
methodology that is better than the nonlinear counterparts available in the literature, in the 
QO . sense of simplicity of implementation and analysis. A scalar control input that excites only 

' one of the states is proposed, and sufficient conditions for the controller gain to stabilize the 

C^) , unstable equilibrium points derived. Numerical simulations confirm the theoretical analysis. 
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1 Introduction 



Keywords: chaos; control of chaos; fractional order systems; stability tests. 



Fractional calculus has been known since the early 17th century [1,. It has been extensively 
applied in many fields, with an overwhelming growth of applications during the last three decades. 
Examples abound in physics [2] , engineering [3] , mathematical biology [4] , finance [5] , life sciences 
[5J, and optimal control [7J. This is due to the fact that, in many applications, approaches based 
on fractional derivatives establish far superior models of engineering systems than the approaches 
based on ordinary derivatives [5] . As mentioned in 1 , there is no field that has remained untouched 
by fractional derivatives. 

Historically, the lack of a physical interpretation of fractional derivatives has been acknowledged 
at the first international conference on the fractional calculus in New Haven (USA), in 1974, by 
including it in the list of open problems [S] . The question was not answered, and therefore repeated 
at the subsequent conferences at the University of Strathclyde (UK) in 1984 [IfJ and at the Nihon 
University (Tokyo, Japan) in 1989 [TT]. The round-table discussion at the conference on transform 
methods and special functions in Varna (1996) showed that the problem was still unsolved [12] . 
Since then, the geometric and physical interpretation of fractional derivatives has been studied by 
several authors [T3J[TJ. An interesting physical discussion about initial conditions of fractional 
order systems is reported in |15j . and their role in control theory is addressed in [16] for Caputo 
fractional derivatives, and in |17j for Riemann-Liouville fractional derivatives. Relation between 
fractional integrals and derivatives and fractal geometry, showing that for some complex systems 



*This is a preprint of a paper whose final and definite form will appear in Mechatronics. Paper submitted 
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substitution of integer-order derivatives by fractional ones results in more accurate and superior 
models, has received special attention [T51ll9j . 

Chaos is an interesting phenomenon in nonlinear dynamical systems that has been developed 
and thoroughly studied over the past two decades. The reader interested in applications of chaos in 
medicine and biology, where fractional calculus has initiated its success and activity in engineering 
applications, is referred to |20 [ l21j. A chaotic system is a nonlinear deterministic system that 
displays complex, noisy-like and unpredictable behavior. The sensitive dependence on the initial 
conditions, and the system parameter variation, is a prominent characteristic of chaotic behavior. 
Here we consider a fractional order chaotic system. The corresponding integer-order system has a 
chaotic behavior for a wide range of parameters. Such integer-order dynamical system is a reduced 
model for a physical system. When such system is implemented using physical electronic devices, 
the environment effects (aging of the elements, temperature, inaccurate values, and so on) on the 
elements appear as a different behavior so that the response predicted by the model does not 
resemble the actual system. Indeed, for such chaotic systems the super sensitivity of the values of 
the elements to tiny changes, cannot be considered in the dynamical equations. The main reason 
is related to the fractal or holed basin of the invariant set of the system. In other words, some 
trace of the system trajectories are seen in the observed coordinate (phase plane) and some of 
them lie in the unseen region that cannot be handled by the classical nonlinear ODE model. For 
this reason, the fractional order system is more appropriate than the classical one [9T fT3llT8"llT9"] . 
This happens to be a frequent situation: for most physical systems which exhibit chaotic behavior, 
the invariant set is not an integer-order dimensional object, and the basin of the trajectories in 
the phase space is a strange attractor field whose Lyapunov dimensions are non-integer. The 
fractional order operators allow to describe the observed behavior by an appropriate kernel in the 
integral. This kernel can be treated as a weighting factor that generates a new response using 
existing vector fields. 

The study of fractional order dynamical systems has attracted an increasing attention in recent 
years due to their great promise as a valuable tool in the modeling of many phenomena 22 . As 
a matter of fact, real world processes generally or most likely are fractional order systems [23] . 
It has been found that fractional order systems possess memory and display more sophisticated 
dynamics when compared to their integer order counterparts, which is of great significance in 
secure communications. On the other hand, due to their potential applications in laser physics, 
chemical reactors, secure communications and economics, a new direction of chaos research has 
emerged in the past two decades to address the more challenging problem of chaos synchronization 
and control [24 271. Recent papers study the control and synchronization of chaotic systems in 
both integer and fractional order cases, applying various control methodologies: nonlinear control 
28,29 , adaptive control [501I5T] . robust control [52"] . fuzzy control [33], and active control [31]. The 
main problem in applying such control methodologies for taming chaos is often their complexity in 
implementation. In this paper we propose a simple linear control mechanism, based on the well- 
known Routh-Hurwitz stability criterion, for a fractional order dynamical system. The system 
shows very rich nonlinear dynamics, including chaos and period doubling bifurcations. A controller 
is proposed, and both analysis and design are studied. We show that a single input can control 
the very complex system. This makes our results simpler than the conventional nonlinear methods 
available in the literature and more feasible to implement. 

The rest of the manuscript is organized as follows. Section [2] briefly presents the necessary 
fractional calculus background. In Section [3] we recall the stability criteria that are used in our 
subsequent analysis and design. Description of an interesting system is presented in Section |U 
while the new results are given in Section [5] where we illustrate our control methodology for 
taming the fractional order chaotic system corresponding to the one of Section [4] Conclusions and 
future directions of research are given in Section [6] 
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2 Fractional Calculus Background 



In this section some necessary mathematical background is presented. For more details see the 
books 



Definition 1 (see, e.g., [10] )• The (left) fractional integral of x of order q, q £ R + , is defined by 

a Dpx{t) ;=J-j\t- S f- X X {s)ds, 



poo 

where T(q) = / e~ z z q ~ 1 dz is the Gamma function. 
Jo 



Definition 2 (see, e.g., The (left) fractional derivative of x of order q, q € R + , in the sense 

of Ricmann-Liouville, is defined by 

^D q x{t) l= D m a D^ m -^x(t) = f ^-^^ j a (t-sr-"- 1 x(s)ds, 

where to € Z + is such that m — 1 < q < m. 

Theorem 3 (see, e.g., |41j). For the fractional Riemann-Liouville derivative and integral one has: 

1. C{ a D- q x{t)} = S -*X{s); 

2. lim Dt q x(t) = D- m x(t), q > 0, m £ 1+ ; 

q—¥m 

3. * L D q c 



"i 1 - — r(l-g) -' 

4. $ L D q tQ Dpx{t) = x(t), q £ R+: 

m—1 

5. C {§ L D q x(t)} = s q X(s) - £ s fc ■ § L D q ~ x(Q), m - 1< q < m, m £ Z+; 

fc=0 

where C denotes the Laplace transform, c a constant, and D~ m the m-folded integral. 

We note that the Ricmann-Liouville differentiation of a constant is not zero (item [3] of The- 
orem [3]) ; while its Laplace transform needs fractional derivatives of the function in initial time 
(item [5] of Theorem [3]). To overcome these "imperfections", the Caputo fractional derivative has 
been introduced. 

Definition 4 (see, e.g., [29] ). The (left) fractional derivative of x of order q, q £ R + , in the sense 
of Caputo, is defined by 

° Dlx® := * L D-^-<»D™x{t) = f^^j j[ (* ~ s) m - q - 1 x^Hs)ds, 

where m £ Z + is such that m — 1 < q < to, D m is the standard differential operator of order m. 
Theorem 5 (see, e.g., |30|). For the fractional Caputo derivative one has: 

1. § D ?c = 0; 

2. §D q D^ q x(t) = § L D q t0 D- q x(t) = x(t), < q < 1/ 

771—1 

5. £{£\D?a;(f)} = s 9 X(s) - £ s^" V^O), to - 1 < q < to, to £ Z+ ; 

k=0 

where c denotes a constant and C the Laplace transform. 

In the next section we study some stability tests for fractional order systems. 
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3 Stability Criteria 

A fractional order dynamical system is usually described by 

\$D q t x{t) = f(x(t),t), m- 1< q < m, m G Z+, t > 0, 
\ $ L D*x(t)] |t=o = x k , k = 0,...,m-l, 

where x(t) e R™ is the vector state at time t, f : R" x K — >• R™ the nonlinear vector field, and 
q = (qi,...,q n ) T the differentiation order vector. If qi = 52 = ■ ■ ■ = q n =: a, we call ([1]) a 
commensurate fractional order dynamical system; otherwise, ([T]) is said to be incommensurate. 

n 

The sum of the orders of all involved derivatives in Eq. ([1]), i.e., ^ qi, is called the effective 

i=l 

dimension of Eq. ((T|) [42] ■ The size of vector x(t) in state space form ([1]), i.e., n, is called the inner 
dimension of system ((T|) [4"3"] . 

Theorem 6 (see |33|). The commensurate order system 

%D?x(t)=Ax(t), x(0)=x o , (2) 

with < a < 1, x(t) € R n , and A € R" xn , is asymptotically stable if, and only if, |arg(A)| > 
for all eigenvalues X of A. Moreover, system ([2J is stable if, and only if, |arg(A)| > for 
all eigenvalues X of A, with those critical eigenvalues satisfying |arg(A)| = having geometric 
multiplicity of one. 

Theorem 7 (see [44]). Consider the following linear fractional order system: 

%D q t x{t)=Ax(t), x(0)=x , (3) 

where x(t) g W 1 , A G R nxn , and q = (q u . . . , q n ) T , < q t < 1, = and gcd(nj,dj) = 1. 
i = 1, . . . , n. If M is the least common multiple of the denominators di, i = 1, . . . , n, then the zero 
solution of ([3]) is globally asymptotically stable in the Lyapunov sense if all roots X of equation 

A(A) = det (diag {X Mq < ) - A) = (4) 

satisfy |arg(A)| > jm- 

Theorem 8 (see [45 ). Let Q = (x*,y*,z*) be an equilibrium solution of |T]) when n — 3 and 
< 51 = q-2 = 53 =: a < 1; and the eigenvalues of the equilibrium point Q for the Jacobian matrix 
^ := dx, \q ^ e 9^ ven by the polynomial A(A) = A 3 + a\X 2 + a^A + 03 = with discriminant 

D(A) = 18aia 2 a 3 + (aia 2 ) 2 - 4a 3 (ai) 3 - 4(a 2 ) 3 - 27(a 3 ) 2 . (5) 

The following holds: 

(i) If D(A) > 0, then a necessary and sufficient condition for the equilibrium point Q to be locally 

asymptotically stable is ai > 0, a 3 > 7 a!a 2 — a 3 > 0. 

(ii) If D(A) < and a\ > 0, a 2 > 0, 03 > 0, then Q is locally asymptotically stable for a < |. 

However, if D(A) < 0, a! < 0, a 2 < 0, a > | , t/ien o/Z roois of Eq. ([5]) satisfy the condition 
|arg(A)|<af. 

(iii) If D(A) < 0, a\ > 0, a 2 > 0, aia 2 — 03 = 0, then Q is locally asymptotically stable for all 
a € (0,1). 

(iv) A necessary condition for the equilibrium point Q to be locally asymptotically stable is 03 > 0. 

(v) If the conditions D(A) < 0, a\ > 0, a 2 > 0, aia 2 — 03 = are satisfied, then the equilibrium 

point Q is not locally asymptotically stable for a = 1 . 
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4 An Interesting System 



In [IB] the following three-dimensional smooth system is proposed and investigated: 

y — ax + byz\ 

cy — xz + z , (6) 
dxy — hz I 

where [x(t), y(t), z(t)] T 6 R 3 is the state vector, and a, b, c, d, and h are some positive constants. 
The system was shown to be chaotic in a wide parameter range, and to have an interesting complex 
dynamical behavior that varies according with the values of the parameters a, 6, c, d, and h. The 
very rich nonlinear dynamics include chaos and period doubling bifurcations. In particular, the 
system generates a two-scroll chaotic attractor for (a, 6, c, d, h) = (3, 2.7, 4.7, 2, 9). 

Chaos may be seen in many real-life nonlinear systems. About ten years ago, several experi- 
mental and theoretical studies have been done to depict the chaotic behavior in various electronic 
systems: nonlinear circuits [47], secure communications [48j[49] , lasers [50], and Colpitts oscil- 
lators [51] . The system ([6j is particularly relevant in mechatronics, where it can be used as a 
carrier producer. Indeed, because of its wide range chaoticity, one of its important applications 
is in secure communication systems. In such a system, a chaotic carrier is used to transmit the 
message signal over a channel. The main motivation for employing such carriers is: (i) complexity 
of the carrier, which increases the security of the modulated signal; (ii) inherent orthonormality 
of the chaotic signals, avoiding the necessity to use in-phase local oscillators as often done in 
telecommunication systems in order to demodulate the original signal at the receiver; (iii) wide 
band signal, which permits the carrier to transmit a wide band message over a noisy channel. 
Motivated by the interesting behavior of ([6j obtained in [46], our main goal is to investigate the 
chaotic dynamics of the corresponding fractional system. This is done in the next section. 



5 Main Results 

Consider a 3D autonomous fractional system 

CD]x(t) = f(x(t)), (7) 

where q = (<7i, <72, <Z3) T is the fractional order of differentiation, x(t) € R 3 is the state vector, and 
/ : R 3 — > R 3 is the nonlinear vector field. Let Q = (x*, x%, x$) be an equilibrium of the system 
([7]), i.e., let f(Q) — 0. We say that Q is a saddle point for ([7J if the eigenvalues of the Jacobian 
matrix J = j£ evaluated at Q are a and b ± jc, where ab < and c ^ 0. A saddle point Q is 
called a saddle point of index 1 if a > and b < 0, and it is called a saddle point of index 2 if 
a < and b > 0. In chaotic systems of Shil'nikov type, scrolls in a chaotic attractor are generated 
only around the saddle points of index 2. Moreover, saddle points of index 1 are responsible only 
for connecting scrolls [52] . 

5.1 System Description 

We are interested in the particular case of ([7]) that corresponds to the commensurate fractional 
order version of (J6j. For that we substitute the standard/integer order derivatives in ([6j by Caputo 
fractional derivatives of order a £ (0, 1): 

f $D?x\ /y-ax + byz\ 

%D? V = [ cy - xz + z . (8) 
q Dfz J \ dxy-hz ) 

System © can be used to model several mechatronic systems. One possible application is to model 
the nonlinear dynamics of a rotor-bearing system with the purpose of diagnosing malfunctions 
and effectively improve the dynamic characteristics of the rotor system. Chu and Zhang analyzed 
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the bifurcation and chaotic motion of a rub-impact rotor system. They found three different 
routes to chaos with an increasing rotating speed [53 . Later, Chu observed very rich forms 
of periodic and chaotic vibrations through experimental verification. These results are of great 
importance to the fault diagnosis of the rub-impact problem [54 . Ehrich studies the bifurcation 
of a bearing-rotor system, identifying a sub-harmonic vibration phenomenon in a rotor dynamic 
system [55]. Goldman and Muszynska analyze the chaotic behavior of a rub- impact rotor using 
numerical emulation and simple experimental verification. They conclude that the rub can lead 
to higher order harmonics, sub-harmonic fractional frequencies, or to chaotic vibrations |56j . Lin 
et al. analyze the nonlinear behavior of rub-related vibration in rotating machinery: the effects 
of the rotating speed, clearance, damping factors, friction coefficients, and boundary stiffness are 
investigated [ST]- Our system ([5]) can also be regarded as a model for a DC-motor with chaotic 
behavior (self-sustained oscillations according to backlash and dead- zone of the gears) [S3]. A 
schematic diagram of a DC drive and its circuits, with separate excitation and controller with 
hysteresis, can be found in [58, Fig. 1]. The states of such system are the current i a in the motor 
armature circuit; the current if in the excitation coil; and the rotor angular speed u r . This real 
life system can be described by ([8]). after normalization and rescale according to the parameters 
and control signals of the circuit, with the correspondence x O uj r , y O if , and 2f>i . 

To find the equilibria of system ([8]), it is enough to equate the right-hand side of ([8]) to zero: 
y — ax + byz = 0, cy — xz + z = 0, dxy — hz = 0. One concludes that the system has 5 equilibria: 

Qi = (o,o,o), 

fd + Va h + VTTa\ -i + VTTa s 

Q2 — 



2d b \ d + VA / 2b 
'd + VA h f-l- VTT~K\ -1-vTTa" 



2d ' b V d+VA J ' 2b J' (9) 
' d- VA h f-l + VTTa\ -i + vTTt\ 



Qs 



2d b V d - VA / 2b 
' d-yfK h f-l - Vl + A\ -1-vT+T" 



2d ' b V d-VA J 7 26 



where A = d 2 + Achd, A = ^ (d + 2ch + VA) , and T = ^ (d + 2ch - VAJ . The Jacobian 
matrix for ([8J, evaluated in an equilibrium point Q = (x*,y*,z*), is given by 



J 



1 + bz* 
c 

dx* 




Our first result gives a necessary condition on the fractional order of differentiation a, for chaos 
to occur. 

Theorem 9 (Necessary condition for occurrence of a chaotic attractor in the fractional order 
system (|8]l). If the fractional order system (J8J exhibits a chaotic attractor, then 

2 /lira (A) h , , 

a > - arctan J (10) 



7T \Re (A) 

for any eigenvalue A of Qi in i = 1, . . . , 5. 

Proof. Assume that the 3D fractional system (JS| displays a chaotic attractor. For every scroll ex- 
isting in the chaotic attractor, the system has a saddle point of index 2 encircled by its respective 
scroll. Suppose that f2 is the set of equilibrium points of the system surrounded by scrolls. A 
necessary condition for the fractional order system ((SJ) to exhibit a chaotic attractor is instability 
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of the equilibrium points in Q |59j . Otherwise, one of these equilibrium points becomes asymp- 
totically stable and attracts the nearby trajectories. According to (U]), this necessary condition is 
mathematically equivalent to 

^-inin{|arg(A l )|}>0, (11) 

where the A, are the roots of dct (diag (\ Mqi A M<?2 A M<?3 ) - J \ Q ) = for all Q e Q. Condition 
(p~0|) follows immediately from (fTTj) . □ 

The nature of the equilibria ([9]) of ((5J may be determined using the corresponding eigenvalues 
A. The following proposition lists the eigenvalues of each equilibrium, when the parameters are 
selected in agreement with Section |U 

Proposition 10. Consider the fractional system (jHJ) of commensurate order a € (0, 1), when the 
five parameters are selected to be (a, b, c, d, h) — (3, 2.7, 4.7, 2, 9). Then, the eigenvalues A for each 
equilibrium Qi in i — 1, . . . , 5, are given as follows. 

(i) Eigenvalues of Q\: —9, —3 and 4.7. 

(ii) Eigenvalues of Q 2 : -11.0247 and 1.8623 ± j6.6831. 

(iii) Eigenvalues ofQ 3 : -11.7856 and 2.2428 ±j6.8580. 

(iv) Eigenvalues of Q 4 : -10.7669 and 1.7335 ± j6.0024. 

(v) Eigenvalues ofQ 5 : -11.6813 and 2.1906 ± j6.1881. 

Proof. Follows by direct computations. □ 

Corollary 11. For the fractional system ((8|) with (a,b,c,d,h) = (3,2.7,4.7,2,9), the equilibria 
Q2, Q3, Qi and Q5 are saddle points of index 2. 

We conclude that if there are some chaotic attractors for (a, b, c, d, h) — (3, 2.7, 4.7, 2, 9), they 
are located around the equilibria Q2, Q3, Q4, Q5. Examining (fTO)) for these equilibria, we obtain: 
a > 0.8270 for Q 2 ; a > 0.7988 for Q 3 ; a > 0.8210 for Q 4 ; and a > 0.7834 for Q 5 . Therefore, by 
choosing a > 0.8270, we ensure that all the eigenvalues remain in the instability region. 

An efficient method for solving fractional order differential equations is the predictor-corrector 
scheme or, more precisely, the PECE (Predict, Evaluate, Correct, Evaluate) technique that has 
been investigated in [Ml EI]- It represents a generalization of the Adams-Bashforth-Moulton 
algorithm. We use the PECE scheme throughout the paper for numerical simulations. 

In Figures Q] to [SJ the initial conditions were selected as (xq, yo, zq) = (5, —2, 1), and only the 
fractional order of differentiation a changes. When a — > 1, our numerical results are in agreement 
with . The numerical simulation of the chaotic attractor for a —> 1 is depicted in Fig. [TJ 

In Fig. [3] we illustrate the chaotic behavior of (|5]) when a = 0.90, and in Fig. [31 Fig. 21 and 
Fig. [5l the behavior of ([8]) is depicted for a — 0.86, a = 0.80, and a = 0.77, respectively. 

As one can see, in Fig. H] and Fig. [S] chaos has diminished, and the trajectories of the corre- 
sponding fractional system converge to their equilibria. 



5.2 Control of the Fractional Order Chaotic System 

Consider the fractional order chaotic system (jH]) of commensurate order a £ (0,1). In order to 
control the system, i.e., force the trajectories to go to the equilibria, we add a control parameter 
u = (wi, it2, W3) as follows: 

o D t x \ /y-ax + byz\ I uA 

D"y = \ cy - xz + z + u a • (12) 
gD?zJ \ dxy-hz J \u 3 J 
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Figure 1: Chaotic attractor of fractional system © with a = 3, b = 2.7, c = 4.7, d = 2, h = 9, 
and initial conditions (xq, yo, Zq) — (5, —2, 1), when a — > 1. 




Figure 2: Numerical results for the fractional order system (jHJ with a = 3, 6 = 2.7, c = 4.7, d = 2, 
h = 9, and initial conditions (x 0l y , z ) = (5, —2, 1), when the fractional order is a = 0.90. 



S 




Figure 3: Numerical results for the fractional order system ([5]) with a — 3, b — 2.7, c — 4.7, d = 2, 
h = 9, and initial conditions (xq, yo, zo) = (5, —2, 1), when the fractional order is a — 0.86. 




Figure 4: Numerical results for the fractional order system ([8]) with a = 3, b — 2.7, c — 4.7, d = 2, 
h = 9, and initial conditions (x ,y , z ) = (5, —2, 1), when the fractional order is a = 0.80. 
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Figure 5: Numerical results for the fractional order system ([8]) with a — 3, b — 2.7, c = 4.7, d = 2, 
h = 9, and initial conditions (x , yo, z ) = (5, —2, 1), when the fractional order is a = 0.77. 

One of the simplest controllers is the state feedback controller, which has a simple structure and 
is easy to implement. Let us take the structure of the controller as a state feedback, as follows: 

= —k\{x — X*), 

-- -k 2 {y- y ), 

-- -k 3 (z-z*). 
In this way (fl~2|) reduces to 



^Dfx\ /y~ax + byz - k\ {x - x*)\ 

CDfy = cy-xz + z- k 2 (y - y*) . (13) 
^Dfz) \ dxy - hz - k 3 (z - z*) ) 

Assume that we want to stabilize one of the equilibria, e.g., Q 2 (using a similar approach, we 
can easily design a control law to stabilize the other unstable equilibria). Next theorem shows 
that by choosing appropriate values for gain ki, we can control the fractional order system (jl3[) . 



Theorem 12. The trajectories of the fractional controlled system (|13j) with a — 3, b — 2.7, c = 4.7, 
d = 2, and h = 9, are driven to the unstable equilibrium point Q 2 = (5.1260, 2.0794, 2.3687) for 
all a G (0, 1), ifk 2 = k 3 =0 and -7.30 < h < 26.53. 

Proof. Computing the Jacobian matrix in the equilibrium point Q = (x*,y*, z*), we obtain 

(-a -ki 1 + bz* by* 
J = -z* c-k 2 -x* + 1 
\ dy* dx* —h — k 3j 

Constituting the characteristic equation A(A) by 

/A + a + 

A(A) = det(AJ - J) = det z* \-c + k 2 x* - 1 1=0, 
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we have 

A(A) = (A + a + ki) [(A - c + k 2 ) (A + h + k 3 ) + dx* (x* - 1)] 

+ (1 + bz*) [z* (X + h + k 3 ) + dy*(x* - 1)] - by* [(-dx*z*) + dy*{X~c + fc 2 )] 
= A 3 + (a + ki + k 2 - c + h + k 3 ) A 2 

+ (dx* (x* ~ 1) + (a + fci)(fe 2 -c + h + k 3 ) + (k 2 - c)(h + k 3 ) + (1 + bz*)z* - bdy* 2 ^j A 

+ (a + h)(dx* (x* - 1) + (fc 2 -c)(h + k 3 )) + (1 + 6z*) (dy*(a;* -l) + z*(h + k 3 )) 
— by* (—dx* z* + dy*(—c + k 2 )) 
= 0. 

(14) 

Based on Theorem [8j if we choose k\, k 2 and k 3 such that 13(A) < 0, a\ > 0, a 2 > and 
aia 2 — a 3 = 0, then Q — (x* , y* . z*) is locally asymptotically stable for all a £ (0, 1), where 

.D(A) = 18 (a + ki + k 2 - c + h + k 3 ) 

x (dx* (x* - 1) + (a + &!)(*& - c + h + ha) + {ha - c)(h + k 3 ) + (1 + 6z*)z* - My* 2 ) 

x ^(a + fci)(tfcc* (s* - 1) + (fc 2 - c)(/i + fc 3 )) + (1 + &O (rfy* - 1) + z*(h + k 3 )) 

- by* {-dx*z* + dy*{~c + fc 2 ))^J 

+ ^(a + ki + k 2 - c + h + k 3 ) (dx* (x* - 1) + (a + h)(k 2 -c+h+ k 3 ) + (k 2 - c)(h + k 3 ) 

i 

+ (1 + bz*)z* - bdy* 2 

4( (a + h)(dx* (x* - 1) + [k 2 - c){h + k 3 )) + (1 + bz*) (dy*{x* - 1) + z*(/i + fe 3 )) 

by*(—dx*z* + dy*(—c + k 2 ))\ (a + ki + k 2 — c + h + k 3 ) 3 
4l dx* {x* - 1) + (a + ki) {k 2 -c + h + k 3 ) + (k 2 -c){h + k 3 ) + [l + bz*)z* - bdy* 2 
27 + fci) (cfe* (x* - 1) + (fc 2 - c)(/i + fc 3 )) + (1 + bz*) {dy*(x* -l) + z*(h + k 3 )) 

3 



- by* (~dx*z* + dy*(~c + k 2 )) 

and the <ij, z = 1,2,3, are found from the coefficients of (fT4"|) . For the parameters that generate 
the 2-scroll attractors, i.e., a = 3, b = 2.7, c = 4.7, = 2, ft = 9, we have Q 2 — (x*,y*,z*) = 
(5.1260, 2.0794, 2.3687). Thus, (TJJ} reduces to 

A(A) = A 3 + (7.3 + h + k 2 + k 3 ) A 2 + (kih 2 + k 2 k 3 + k x k 3 + 21.8733) A + 530.6404 

- 0.0002fc! + 3.6509fc 2 + 3.4177/c 3 + 3fc 2 fc 3 + 9fc!fc 2 - 4.7feife 3 + hk^. 

If we choose k 2 = k 3 = 0, then 

A(A) = A 3 + (7.3 + fti) A 2 + 21.8733A + 530.6404 - 0.0002/ci. 
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Figure 6: Numerical results for the controlled system (fT3f with a — 3, b — 2.7, c = 4.7, d = 2, 
ft, = 9, fci = 16.96, k-2 = &3 = 0, when the fractional order is a — 0.90 and initial conditions are 
(x Q ,y Q ,z ) = (5,2,2). 

For this characteristic polynomial we have: 

D(A) = 18 (7.3 + fci) (21.8733) (530.6404 - 0.0002fci) + ((7.3 + fci) 21.8733) 2 

- 4 (530.6404 - 0.0002fci) (7.3 + fci) 3 - 4(21. 8733) 3 - 27(530.6404 - 0.0002fci) 2 . 

The conditions D(A) < 0, ai > 0, <22 > and 0102—03 = arc satisfied for k\ € (—7.30, 26.53). □ 

In Fig. [6] we have chosen fci = 16.96 and the initial conditions to be (xo, yo, zq) — (5, 2, 2). As 
can be seen from the figure, all states converge to their equilibria. 

In order to get a faster response, one can easily consider the other two gains (fc2 and k^) in the 
control law. However, this introduces difficulties in the real implementation of the control system. 
We also note that in designing the stabilizing controllers, because we are utilizing the linearized 
version of the nonlinear system around its equilibrium, we should select the initial conditions near 
the corresponding equilibrium. 

6 Conclusions and Future Work 

Chaotic fractional order systems have an inherent potential in mechatronic applications, particu- 
larly in secure telecommunication systems where the main part of a transmitter-receiver config- 
uration is the synchronization between master and slave blocks. In this article we analyzed the 
dynamical behavior of a novel fractional order chaotic system. The chaotic system generalizes 
the recent integer order system introduced in [3B]. The local stability of the equilibria, using 
the fractional Routh-Hurwitz conditions, was studied. Furthermore, using Matignon's stability 
criteria, the system was shown to be chaotic in a wide parameter range, and to have an inter- 
esting complex dynamical behavior that varies according with the values of the parameters a, 6, 
c, d, and h. The very rich nonlinear dynamics include chaos and period doubling bifurcations. 
Moreover, we derived a lower bound of the fractional order of differentiation for the system to 
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remain chaotic. Analytical conditions for linear feedback control have been achieved. Our anal- 
ysis is valid in a general format, in which all gains k\, k%, and k% are considered. However, and 
despite the complexity and the wide range of varieties, it has been shown that the system can be 
controlled by a single state-feedback controller. This possibility of stabilizing the system locally 
using only controller k\ provides a simple and easy way to control the chaos, which can be crucial 
in a real implementation. Our state feedback approach can be applied to various chaotic fractional 
order systems. In particular, we claim that the techniques here developed can be used, without 
fundamental changes, in the synchronization of two fractional order chaotic systems in a secure 
telecommunication system. This is under study and will be addressed elsewhere. 

Our state feedback approach for controlling the system is based on a Routh-Hurwitz analysis 
and cannot treat constraints on the actuating signal. To consider such saturation constraints 
on the control signal, a performance index should be defined. Minimizing the index subject to 
some constraints on the state and control signals cannot be considered analytically for the method 
here implemented. A direction for future research is to investigate how one can obtain a global 
stabilizing controller. To the best of our knowledge, the global stability of fractional systems is an 
interesting open question and available results reduce to those of [62l[63]. For different stability 
concepts than the Lyapunov one adopted in our work, we refer the reader to [64j . 
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